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Abstract 

We  study  the  overall  thermal  conductivity  of  composites  under 
three  different  heating  scenarios:  (i)  a  laser  pulse  heat  source,  (ii) 
a  preheated  composite  sample,  and  (iii)  a  continuous  heat  source. 


1  Introduction 

Adhesives  such  as  epoxies,  gels,  and  greases  have  numerous  commercial  and 
industrial  applications,  such  as  computers,  machinery,  and  home  appliances. 
These  adhesives  are  typically  very  poor  conductors.  Hence,  filler  particles, 
with  higher  thermal  conductivities,  such  as  diamond  dust,  carbon  fibers,  or 
aluminum  particles,  are  added  to  create  a  composite  material  that  is  a  better 
thermal  conductor  than  the  original  material.  The  filler  particles  introduce 
several  design  considerations,  including  the  choice  of  particle,  the  particle 
geometry,  and  the  size  and  shape  of  the  particles. 

In  [1]  the  authors  present  a  mathematical  model  to  describe  the  two 
dimensional  heat  transfer  through  a  composite  silicone.  The  physical  model 
involves  a  thin  slice  of  composite  material  with  a  pulsed  heat  source  (laser) 
at  one  end,  a  heat  sink  at  the  other  end  (ambient  air),  and  no  heat  loss 
through  the  sides.  See  Figure  1  for  a  depiction  of  the  physical  model.  The 
corresponding  mathematical  model  is  based  on  the  transient  heat  equation 
with  insulated  boundary  conditions  on  the  sides,  a  flux  due  to  the  heat  source 
at  one  end,  and  Newton  cooling  at  the  other  end.  Let  u(t,z)  represent  the 
composite  silicone  temperature  at  a  given  time  t  and  coordinate  z.  Let  k(z ), 
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Figure  1:  Two  dimensional  heat  transfer  model 


p(z),  and  cp(z)  represent  the  thermal  conductivity,  density,  and  specific 
heat,  respectively,  of  the  composite  at  coordinate  z.  Then  the  following 
system  describes  the  temperature  in  the  composite  sample  for  zEfl  — 
[-f,f]x[-f,f]andtG[0,T], 


p(z)cp(z)  —  (t,z)  =  V  •  (k(z)Vu(t,z)),  u(0,z)  =  u0(z),  zEVt 

Hz)^(^z)\n  =  o 

du 

k{Z)  (^5  Z)  1 72  —  MToo  —  U{tt  X))  I72  (1) 

k(Z)^(^Z)\l3  =  0 

*(*)§£(M)It.  =  so{t)- 


Here  denotes  the  ambient  temperature,  h  is  the  Newton  cooling  constant 
(or  heat  transfer  coefficient),  So(t)  is  the  flux  due  to  the  heat  source,  and 
uo  is  the  initial  temperature  distribution  in  the  composite  silicone.  The 
boundaries  are  described  by  the  parameterized  7 ,,  where 


7i(*)  =  (y,*):*e[- 1,|], 
73  W  = 


72(5)  =  (i,|):*e  [-y, 

74(s)  =  (a,-y)  :  a  e  by* 
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and  c?0  =  71U72U73U74.  See  Figure  2  for  a  depiction  of  the  two  dimensional 
slice  with  indicated  boundaries.  Let  0.5  denote  the  region  occupied  by  the 
silicone  and  Op  denote  the  region  occupied  by  the  filler  particles.  Then  0  = 
0.5  U  Op  and  0  =  0.,  n  Op.  Note  that  the  coefficients  of  the  partial  differential 
equation,  k(z ),  p(z),  and  cp(z),  are  piecewise  constant.  For  example,  k(z )  = 
ks  for  z  E  fls  and  k(z )  =  kp  for  z  E  Op. 
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Figure  2:  2D  composite  adhesive  slice 


In  [1]  the  authors  show  this  model  is  well-posed,  i.e.,  there  exist  unique 
weak  solutions  which  are  continuously  dependent  on  the  initial  conditions, 
forcing  function  and  parameters.  They  also  develop  mathematical  models 
that  describe  the  sensitivity  of  temperature  on  the  silicone  and  filler  particle 
thermal  conductivities  and  establish  the  well-posedness  of  the  sensitivity 
models.  Their  computational  investigations  involve  a  uniform  arrangement 
of  circular  filler  particles  within  the  silicone.  From  these  investigations, 
the  authors  conclude  that  the  model  is  not  very  sensitive  to  filler  particle 
conductivity. 

In  this  note  we  continue  the  computational  investigations  of  [1]  to  further 
understand  the  heat  transfer  through  a  composite  silicone.  We  concentrate 
our  presentation  on  a  composite  material  of  fixed  size  and  consider  the  effect 
of  particle  size  and  arrangement.  Our  investigations  explore  three  experi¬ 
mental  scenarios.  The  first  scenario  involves  a  laser  pulse  heat  source  as 
described  in  [1]  and  above.  This  scenario  imitates  a  heat  surge  above  a 
steady  state  condition.  Solutions  for  this  scenario  are  highly  transient,  par¬ 
ticularly  with  respect  to  the  sensitivity  analysis.  Hence,  for  a  second  scenario 
we  consider  an  initially  heated  adhesive  with  three  insulated  boundaries  and 
the  fourth  boundary  exposed  to  cooling  (ambient  air).  The  third  scenario 
involves  a  continuous  heat  source  and  focuses  on  steady  state  solutions.  This 
scenario  imitates  a  situation  wherein  silicone  in  a  computer  motherboard  is 
in  the  presence  of  continuous  power  source  such  as  during  regular  computer 
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usage. 

The  numerical  experiments  are  outlined  below.  In  Section  2  we  report  on 
experiments  with  a  laser  pulse  heat  source.  Much  of  the  detail  and  discussion 
is  included  in  this  first  to  better  understand  the  basic  heat  transfer  behavior. 
In  Sections  3  and  Section  4  we  consider  the  preheated  and  continuous  heat 
source  scenarios,  respectively.  The  presentations  in  these  later  sections  are 
similar  to  Section  2  and  their  results  are  intended  for  comparison.  Some 
conclusions  are  summarized  in  Section  5. 

Our  numerical  results  will  show  that,  in  general,  a  smaller  particle  size 
and  larger  volume  fraction  can  slightly  improve  the  heat  transfer  through 
composite  adhesives.  However,  larger  gains  may  be  realized  by  increasing 
the  silicone  thermal  conductivity  in  the  case  of  a  continuous  or  pulsed  heat 
source. 


2  Pulse  Heat  Source  Experiments 


We  consider  the  basic  model  presented  above,  i.e.,  insulated  boundaries  on 
the  sides,  a  laser  pulse  heat  source  on  one  end,  and  Newton  cooling  on  the 
other  end.  The  size  of  the  adhesive  sample  is  107.67  x  546,  i.e.,  ci  =  107.67 // 
and  C2  =  546//.  The  laser  pulse  heat  source  is  a  flux  term  of  the  form 


So 


Si 

0 


o  <  t  <  tp 

tp  <c  t, 


where  tp  —  0.000330  seconds.  The  model  is  numerically  solved  on  a  tri¬ 
angular  mesh  generated  by  the  Matlab  PDE  Toolbox  using  the  Delaunay 
triangulation  algorithm.  We  integrate  system  (1)  on  the  interval  t  G  [0,  3.46] 
seconds  using  a  finite  element  method  from  Matlab’s  PDE  Toolbox.  See  [1] 
for  details  about  the  source  flux  or  the  use  of  the  PDE  Toolbox. 

In  working  with  transient  solutions  (i.e.,  a  pulse  source)  the  temper¬ 
ature  profiles,  flux,  and  in  particular,  sensitivity  results,  are  highly  time 
dependent.  Additionally,  the  time  frame  for  rapid  changes  is  affected  by  the 
silicone  thermal  conductivity.  To  better  understand  the  basic  heat  transfer 
behavior  we  first  look  at  an  adhesive  sample  without  particles.  The  adhe¬ 
sive  will  have  the  material  characteristics  of  either  silicone  or  filler  particle. 
We  then  look  at  filled  composite  adhesives,  with  either  a  uniform  or  ran¬ 
dom  particle  arrangement.  Finally,  we  consider  the  effect  of  particle  size, 
placement  and  volume  fraction  on  heat  transfer. 
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2.1  Without  Particles 


To  investigate  the  effect  of  changing  the  thermal  conductivity,  we  con¬ 
sider  an  unfilled  adhesive  sample.  We  work  with  two  cases:  a  sample 
with  thermal  conductivities  in  the  silicone  range  (k  —  0.02  —  1.02 W/mK), 
and  a  sample  with  thermal  conductivities  in  the  filler  particle  range  (k  — 
100  —  1100 W/mK).  The  density  and  specific  heat  parameters  are  set  ac¬ 
cording  to  the  material  in  each  case.  The  model  parameters  are  listed  in 
Table  1. 


UW/mK) 

dg/cm6 ) 

cp/gK) 

Too{K) 

h 

S(W/m') 

0.02  -  1.02 

1 9 

1.55 

296.15 

350 

4.32aT0y 

100  -  1100 

2.7 

0.90 

296.15 

350 

4.32aT0y 

Table  1:  Model  parameters 

In  Figure  3  we  plot  the  temperature  in  the  sample  at  various  times 
during  the  time  integration  for  two  different  silicone  thermal  conductivi¬ 
ties.  The  first  column  of  plots  indicate  temperatures  for  a  sample  with 
k  —  0.02 W/mK]  the  second  column  indicates  temperatures  for  a  sample 
with  k  —  1.02 W/mK.  At  t  —  0.0115  seconds,  the  pulse  may  be  clearly 
identified  in  one  end  of  the  slice.  After  t  —  1.1418  seconds,  the  pulse 
is  distributed  through  the  slice,  and  at  the  end  of  the  time  integration 
( t  —  3.46s),  most  of  the  heat  is  convected  away  by  the  Newton  cooling 
boundary  condition.  Note  the  difference  in  temperature  range  for  a  slice 
with  low  (0.02 W/mK)  versus  high  (Uy2W/rnK)  silicone  thermal  conduc¬ 
tivity.  Increasing  the  silicone  thermal  conductivity  results  in  a  large  im¬ 
provement  in  heat  transfer. 

We  next  consider  a  sample  with  thermal  conductivities  in  the  range  of 
filler  particles.  In  Figure  4  we  plot  the  temperature  in  the  sample  at  var¬ 
ious  times  for  two  different  thermal  conductivities.  The  first  column  plots 
temperatures  with  k  —  100 W/mK]  the  second  column  plots  temperatures 
with  k  —  1100 W/mK.  In  each  case,  after  only  t  —  0.0115  seconds,  much  of 
the  pulse  has  moved  through  the  sample.  There  is  relatively  little  difference 
in  temperature  between  plots  with  low  and  high  particle  thermal  conduc¬ 
tivity.  Hence,  in  this  experiment,  changing  the  particle  conductivity  alone 
does  not  greatly  affect  the  resulting  heat  transfer  within  the  sample.  Note 
the  difference  in  temperature  ranges  for  a  particle  sample  (Figure  4)  versus 
a  silicone  sample  (Figure  3). 

Next  we  compute  the  average  temperature  in  a  slice  at  a  particular 
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Figure  4:  Filler  particle  sample 
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Average  temperatures:  silicone 


Figure  5:  Average  temperatures  at  t  —  1.14s  and  t  —  3.46s. 


time  and  compare  average  temperatures  over  a  range  of  adhesive  thermal 
conductivities.  We  compute  the  average  temperature  as  follows 

|A|_1  /  u(t,z-,k)dz, 

Ja 

where  A  represents  the  area  of  the  adhesive.  In  Figure  5  we  plot  the  average 
temperature  in  the  slice  for  k  in  the  range  of  silicone  conductivities,  and  for  k 
in  the  range  of  filler  particle  conductivities.  The  solid  line  represents  average 
temperatures  at  t  —  1.1418  seconds  into  the  time  integration;  the  dashed  line 
represents  average  temperatures  at  the  final  time  t  —  3.46  seconds.  Recall 
that  the  ambient  temperature  is  =  296.15A4  Clearly,  increasing  silicone 
conductivity  is  more  effective  in  lowering  the  average  temperature  in  the 
slice,  though  the  effect  levels  off  at  higher  silicone  conductivities.  Increasing 
particle  conductivity  has  no  visible  effect. 

The  basic  heat  transfer  mechanisms  present  are  conduction  and  convec¬ 
tion.  Heat  is  transfered  through  the  slice  by  conduction;  heat  is  transfered 
from  the  slice  to  the  ambient  air  by  convection  due  to  the  Newton  cool¬ 
ing  boundary  condition  at  72.  From  the  boundary  condition,  the  flux  at 
72  is  proportional  to  the  temperature  gradient  at  72.  In  Figure  6  we  plot 
the  flux  at  72,  averaged  over  72,  as  a  function  of  time  for  various  thermal 
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Figure  6:  Average  flux  at  72. 


conductivities.  Specifically,  we  plot 

77  7  du  .  T72  .1  1-1  f  du ,  \  , 

=-k\l2\-  J^(t,S)ds. 

The  top  plot  gives  the  average  flux  for  thermal  conductivities  in  the  sil¬ 
icone  range,  k  —  ks  £  [0.02, 1.02] W/mK]  the  bottom  plot  displays  the 
average  flux  for  thermal  conductivities  in  the  filler  particle  range,  k  —  kp  G 
[100, 1100]  W/mK.  For  conductivities  in  the  particle  range,  increasing  the 
conductivity  has  little  effect  on  the  average  flux  at  72.  Conversely,  for  con¬ 
ductivities  in  the  silicone  range,  increasing  the  conductivity  does  increase 
the  average  flux  at  72.  This  increase  saturates  with  the  larger  silicone  con¬ 
ductivities  due  to  the  Newton  cooling  boundary  condition  at  72  and  the 
cooling  rate  (h). 
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The  constant  h  depends  on  several  factors:  the  adhesive  properties  ( k ,  /;, 
p1  Cp).  the  system  geometry,  the  flow  velocity,  the  value  of  the  characteristic 
temperature  difference  at  the  boundary,  as  well  as  surface  temperature  dis¬ 
tribution.  For  our  experiments,  the  value  of  h  mimics  ambient  air  cooling. 
The  use  of  fans  or  other  devices  could  result  in  a  faster  or  slower  cooling 
rate,  respectively. 

We  now  look  at  results  from  the  sensitivity  equations.  Let  w{t,z)  — 
||  (t.  z )  represent  the  sensitivity  of  temperature  to  thermal  conductivity.  See 
[1]  for  the  sensitivity  equations  and  their  derivation.  Since  the  temperature 
solutions  are  transient,  the  sensitivities  are  highly  time-dependent.  In  [1] 
the  authors  plot  the  sensitivity  at  72  (average  sensitivity  or  relative  average 
sensitivity)  as  a  function  of  conductivity  at  various  instances  in  time.  Instead 
we  plot  sensitivity  at  72,  averaged  over  72,  as  a  function  of  time  for  various 
adhesive  thermal  conductivities.  More  specifically,  we  compute  the  average 
sensitivity 

W2{t]k)  —  |72|_1  /  w(t,s]k)ds 

J  72 

for  k  in  the  silicone  range  or  the  filler  particle  range.  The  top  plot  in  Fig¬ 
ure  7  displays  the  average  sensitivity  to  silicone  thermal  conductivity  at  72; 
the  bottom  plot  displays  the  average  sensitivity  to  particle  thermal  con¬ 
ductivity  at  72.  A  positive/negative  average  sensitivity  indicates  that  the 
temperature  is  increasing/decreasing  (on  average)  with  respect  to  k  at  72. 
Recall  that  the  slice  is  initially  at  ambient  temperature.  As  the  pulse  dis¬ 
tributes  through  the  slice  the  temperature  at  the  sink  boundary  72  initially 
increases.  Increasing  the  thermal  conductivity  has  the  effect  of  transferring 
more  heat  from  the  source  end  to  the  sink  end,  some  of  which  is  not  imme¬ 
diately  convected  away  and  results  in  a  higher  temperature  at  72.  As  heat  is 
convected  from  the  slice,  temperatures  eventually  decrease  at  72.  Increasing 
the  thermal  conductivity  results  in  more  effective  heat  transfer  to  72  and 
lower  temperatures.  For  conductivities  in  the  particle  range,  the  decrease 
in  temperature  occurs  almost  immediately;  for  conductivities  in  the  silicone 
range,  temperatures  initially  increase  and  then  decrease. 

Note  the  different  ranges  for  average  sensitivity  to  silicone  conductivity 
versus  particle  conductivity,  indicating  that  heat  transfer  is  more  sensitive  to 
changes  in  silicone  conductivity  than  particle  conductivity.  Also  note  that 
for  t  >  Is,  temperature  is  insensitive  to  both  particle  or  silicone  thermal 
conductivity  (except  for  the  lowest  silicone  conductivity) . 
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Test  case:  all  silicone,  v^(t;k  ) 


Figure  7:  Average  sensitivity  at  boundary  72. 

2.2  With  Particles 

In  this  section  we  include  filler  particles  within  the  adhesive.  The  particles 
are  circular  and  uniform  in  size  with  radius  12.07/r;  the  particles  are  alu¬ 
minum  with  fixed  thermal  conductivity  kp  —  ‘1Y1W jinK.  Since  the  silicone 
polymer  is  assumed  to  wet  well,  we  assume  a  minimum  distance  between 
particles  and  that  no  particles  touch  the  boundary  of  the  composite  silicone 
slice.  The  size  of  the  silicone  slice  is  the  same  as  above. 

We  consider  two  different  particle  arrangements:  a  uniform  arrangement 
of  thirty- two  particles  aligned  in  two  rows  (UNI) ,  and  a  random  arrangement 
of  thirty- two  particles  (RAND).  Both  arrangements  constitute  a  volume 
fraction  of  25%.  The  results  reported  in  [1]  refer  only  to  a  uniform  particle 
arrangement,  however  a  random  arrangement  is  more  realistic.  See  Figure 
8  for  the  particle  arrangements.  The  purpose  of  these  experiments  is  to 
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Figure  8:  Composite  adhesive  particle  arrangement  vf  25% 


determine  the  effect  of  different  particle  arrangements,  and  also  to  compare 
with  results  from  silicone  slices  in  the  previous  subsection. 

In  Figure  9  we  plot  the  temperature  profiles  resulting  from  the  uniform 
and  random  particle  arrangements  at  the  intermediate  time,  t  —  1.1418s. 
The  top  two  plots  provide  the  profiles  for  the  uniform  and  random  arrange¬ 
ments  with  a  low  silicone  thermal  conductivity  ( ks  —  Q.(y2W/rnK):  the 
bottom  two  plots  provide  profiles  for  each  arrangement  with  a  high  sili¬ 
cone  thermal  conductivity  ( ks  —  1.02 W/mK).  Again  the  particle  thermal 
conductivity  is  set  for  aluminum  with  kp  —  217 W/mK.  There  are  slight 
differences  in  profiles  due  to  the  particle  arrangements  (i.e. ,  heated  rows 
versus  irregular  warm  spots),  however,  the  overall  temperatures  for  each 
arrangement  appear  similar. 

In  Figure  10  we  plot  the  average  temperature  in  the  uniform  composite 
slice  (UNI),  the  random  composite  slice  (RAND),  and  the  all  silicone  slice 
(SIL)  at  t  —  1.1418s  and  at  the  final  time  t  —  3.46s  for  a  range  of  sili¬ 
cone  thermal  conductivities.  The  average  temperatures  for  the  composite 
slices  (UNI  and  RAND)  are  nearly  identical  at  either  time.  At  the  inter¬ 
mediate  time  ( t  —  1.1418s),  the  composite  silicone  slices  have  lower  average 
temperatures  than  the  all  silicone  slice  across  the  range  of  silicone  thermal 
conductivities.  In  other  words,  the  composite  slices  are  more  effective  in 
transferring  heat  through  and  out  of  the  adhesive  slice.  This  is  also  true  at 
the  final  time  ( t  —  3.46s),  however,  the  average  temperatures  saturate  as 
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Figure  9:  Temperature  profiles  for  composite  adhesives 


Average  temperatures:  t=1 .1 41 8s,  kp=217W/mK 


Figure  10:  Average  temperature  for  composite  adhesives,  vf  25% 


ks  increases  since  most  of  the  heat  has  convected  away.  Little  improvement 
in  heat  transfer  is  realized  by  increasing  the  silicone  thermal  conductivity 
beyond  0.3 W/mK. 

2.3  Fixed  Volume  Fraction 

In  the  previous  subsection  the  uniform  versus  random  particle  placement 
did  not  noticeable  affect  the  resulting  temperatures  from  a  pulse  source.  In 
this  section  we  investigate  in  more  detail  the  effect  of  particle  placement  and 
particle  size  on  heat  transfer. 

We  fix  the  volume  fraction  at  25%  and  vary  the  size  and  placement  of  the 
particles.  Again,  we  work  with  aluminum  filler  particles  ( kp  —  2l7W/mK). 
All  slices  contain  a  uniform  particle  placement:  3  slices  with  particles  ar¬ 
ranged  in  two  rows,  3  slices  with  particles  arranged  in  1  row.  In  each  case  the 
number  of  particles  is  inversely  proportional  to  the  size  of  the  particles.  See 
Figure  11  for  the  different  particle  arrangements  with  25%  volume  fraction. 

To  better  observe  trends  in  heat  transfer,  we  plot  the  temperatures 
through  the  middle  of  each  slice  ( x  —  0)  at  the  time  t  —  1.1418s  for  two 
different  silicone  thermal  conductivities  in  Figure  12.  The  top  two  plots  give 
temperature  profiles  for  a  low  silicone  conductivity  ( ks  —  0.02 W/mK)]  the 
bottom  two  plots  show  temperature  profiles  for  a  high  silicone  conductivity 
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Figure  11:  Uniform  geometries  with  25%  volume  fraction 


Figure  12:  Temperature  profiles  at  x  —  0,  25%  volume  fraction 
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Average  temperatures,  t=1 .1418s,  k=217  W/mK 


Figure  13:  Average  temperatures,  25%  volume  fraction 


(ks  —  1.02 W/mK).  Temperatures  are  plotted  at  an  intermediate  time  since 
most  of  the  laser  pulse  has  dissipated  away  by  the  end  of  the  time  integra¬ 
tion.  The  differences  in  temperature  profiles  for  each  of  the  different  uniform 
geometries  are  slight.  However,  the  heat  transfer  for  the  eight  particle  case 
is  more  efficient  with  particles  aligned  in  one  row  as  opposed  to  two  rows. 
Placement  of  the  particles  in  one  row,  with  the  particles  closer  together,  acts 
as  a  channel  to  direct  the  heat  transfer  through  the  filler  particles. 

In  Figure  13  we  plot  the  average  temperature  for  each  of  the  uniform  par¬ 
ticle  arrangements  at  the  time  t  —  1.1418s  across  a  range  of  silicone  thermal 
conductivities.  For  reference  we  have  also  included  the  average  temperature 
in  the  all  silicone  slice  (SIL)  at  the  same  time  across  the  same  range  of 
conductivities.  There  are  slight  differences  in  the  average  temperatures  be¬ 
tween  the  various  uniform  geometries  of  25%  volume  fraction,  though  they 
all  remain  close  together.  In  general  the  composite  slices  with  smaller,  more 
numerous  particles  produce  lower  average  temperatures.  All  the  composite 
slices  obtain  lower  average  temperatures  than  the  silicone  slice  by  approxi¬ 
mately  2  K. 
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Figure  14:  Random  geometries  with  varied  volume  fraction 

2.4  Varied  Volume  Fraction 

In  this  section  we  look  at  the  effect  of  volume  fraction  on  heat  transfer 
through  a  composite  adhesive.  Here  we  fix  the  particle  size  ( R  —  10.50//) 
and  vary  the  number  of  particles  (volume  fraction)  in  the  composite  adhe¬ 
sive.  The  particles  are  arranged  randomly  to  form  volume  fractions  of  25%, 
35%,  45%,  and  55%.  See  Figure  14  for  the  random  particle  arrangements. 
The  particles  are  aluminum  ( kp  —  217 W/mK).  The  size  of  the  slice  and 
particles  are  smaller  than  in  the  previous  section  to  ease  the  domain  gener¬ 
ation  performed  by  MATLAB  for  the  high  volume  fraction  domains.  Here 
ci  =  80.375//  and  c 2  =  409.5 fi. 

We  are  looking  for  trends  in  heat  transfer  behavior  as  we  increase  volume 
fraction.  From  the  experiments  with  fixed  volume  fraction,  there  appears  to 
be  a  slight  benefit  to  smaller,  closer  particles.  The  expectation  is  that  the 
higher  the  volume  fraction,  the  more  effective  the  heat  transfer. 

In  Figure  15  we  plot  the  temperature  profiles  through  the  middle  of  the 
slice  ( x  —  0)  at  the  time  t  —  1.1418s  for  two  different  silicone  thermal  con¬ 
ductivities.  The  top  plot  displays  temperature  profiles  for  the  low  silicone 
conductivity  ( ks  —  0.02 W/mK)]  the  bottom  plot  displays  temperature  pro¬ 
files  for  the  high  silicone  conductivity  ( ks  —  1 ,02W/mK).  Both  plots  display 
the  same  trend:  larger  volume  fractions  result  in  lower  temperatures.  The 
total  temperature  difference  is  less  than  1 K,  however,  for  slices  with  high 
silicone  thermal  conductivity. 
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Figure  17:  Initially  heated  2D  composite  adhesive  slice 


This  trend  is  more  clearly  seen  in  Figure  16  where  we  plot  the  average 
temperatures  for  each  slice  at  the  time  t  —  1.1418s  across  a  range  of  silicone 
thermal  conductivities.  Again,  we  include  the  average  temperature  for  the 
all  silicone  slice  (SIL)  for  the  same  range  of  silicone  conductivities.  All 
composite  slices  have  a  lower  average  temperature  than  the  all  silicone  slice; 
The  higher  the  volume  fraction,  the  lower  the  average  temperature.  Note, 
however,  that  the  difference  in  average  temperatures  from  25-55%  volume 
fraction  is  only  1.5 K.  Also  note  that  all  average  temperatures  saturate, 
at  ks  >  0.3W/mK  for  composite  slices  and  at  ks  >  0.7 W/mK  for  the 
silicone  slice.  Prior  to  saturation,  however,  an  increase  in  silicone  thermal 
conductivity  results  in  lower  average  temperatures.  Hence,  another  means 
of  improving  heat  transfer  may  be  accomplished  by  increasing  the  silicone 
thermal  conductivity  up  to  the  saturation  point. 

3  Heated  Slice  Experiments 

In  this  section  we  work  with  a  slightly  different  experimental  set-up  to  better 
avoid  the  transient  nature  of  solutions  in  the  previous  section.  Rather  than 
using  a  laser  pulse  heat  source,  the  slice  is  pre-heated  to  a  temperature 
of  650/F.  Three  sides  of  the  sample  are  insulated  (71,  73,  and  74);  heat 
convects  from  the  sample  through  a  Newton  cooling  boundary  condition  on 
the  fourth  side  (72).  See  Figure  17. 
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Figure  18:  Temperature  profiles  at  x  —  0,  25%  volume  fraction 


In  the  following  subsections  we  investigate  the  following  effects  on  heat 
transfer:  varying  particle  size  and  placement  for  a  fixed  volume  fraction, 
and  varying  volume  fraction  for  a  fixed  particle  size. 

3.1  Fixed  Volume  Fraction 

Here  we  investigate  the  effect  of  particle  size  and  placement  on  heat  transfer 
through  a  heated  composite  adhesive.  We  fix  the  volume  fraction  and  vary 
the  size  and  placement  of  the  particles.  We  work  with  the  same  six  uniform 
configurations  of  25%  volume  fraction  as  in  the  previous  section:  3  samples 
with  particles  arranged  in  two  rows,  3  samples  with  particles  arranged  in  1 
row.  See  Figure  11  for  the  different  particle  arrangements.  Again,  the  filler 
particles  are  aluminum  (kp  —  217 W/mK). 

In  Figure  18  we  plot  the  temperature  profiles  through  the  middle  of 
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Average  temperatures:  t=3.46s,  kp=217W/mK 


Figure  19:  Average  temperature  for  heated  slices  with  25%  volume  fraction 


each  slice  (x  —  0)  at  the  final  time  t  —  3.46s  for  two  different  silicone 
thermal  conductivities.  The  top  two  plots  give  temperature  profiles  for  a 
low  silicone  conductivity  ( ks  —  0.02 W/mK)  slice;  the  bottom  two  plots  show 
temperature  profiles  for  a  high  silicone  conductivity  (ks  —  l.02W/mK)  slice. 
Temperatures  are  plotted  at  the  final  time,  rather  than  an  intermediate 
time,  since  there  is  plenty  of  heat  remaining  in  the  composite  slice.  The 
temperature  trends  are  more  pronounced  than  in  the  pulsed  heat  source 
scenario.  In  both  the  one  row  and  two  row  configurations,  the  smaller,  closer 
particle  arrangements  result  in  better  heat  transfer,  i.e.,  low  temperatures. 
Also,  for  the  eight  particle  case,  heat  transfer  is  more  efficient  when  the 
particles  are  aligned  in  one  row  as  opposed  to  two  rows. 

In  Figure  19  we  plot  the  average  temperature  for  each  of  the  uniform 
particle  arrangements  at  the  final  time  t  —  3.46s  across  a  range  of  silicone 
thermal  conductivities.  Also  included  is  the  average  temperature  in  the 
silicone  slice  (SIL)  at  t  —  3.46s  across  the  same  range  of  conductivities. 
There  are  slight  differences  in  the  average  temperatures  between  the  various 
uniform  geometries  of  25%  volume  fraction,  though  the  temperature  ranges 
are  larger  than  in  previous  experiments.  In  general  the  composite  slices 
with  smaller,  more  numerous  particles  produce  lower  average  temperatures. 
The  greatest  difference  in  average  temperatures  is  between  the  composite 
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Figure  20:  Temperature  profiles  at  x  —  0.  varied  volume  fraction  volume 
fraction 

slices  and  the  all  silicone  slice.  With  a  low  silicone  thermal  conductivity,  the 
composite  slices  have  a  lower  average  temperature  than  the  silicone  slice. 
However,  for  ks  >  0.02 W/mK,  the  silicone  slice  produces  a  lower  average 
temperature  at  the  final  time.  This  cross-over  point  depends  on  the  heat 
transfer  coefficient  h. 

The  average  temperature  plot  in  Figure  19  indicates  that,  for  large  sili¬ 
cone  conductivity,  the  silicone  slice  appears  to  transfers  the  heat  through  the 
slice  more  effectively  than  the  composite  slices  with  similar  silicone  thermal 
conductivity.  This  result  does  not  match  results  from  Section  2. 

3.2  Varied  Volume  Fraction 

In  this  section  we  look  at  the  effect  of  volume  fraction  on  heat  transfer 
through  a  heated  composite  adhesive.  We  fix  the  particle  size  ( R  —  10.50//.) 
and  vary  the  number  of  particles  (volume  fraction)  in  the  composite  ad¬ 
hesive.  The  particles  are  aluminum  (kp  —  2YIW/mK)  and  are  arranged 
randomly  to  form  volume  fractions  of  25%,  35%,  45%,  and  55%.  See  Figure 
14  for  the  random  particle  arrangements. 
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Average  temperatures:  t=3.46s,  kp=217W/mK 


Figure  21:  Average  temperatures,  varied  volume  fraction 


Again,  we  look  for  trends  in  heat  transfer  as  we  increase  the  volume 
fraction.  In  Figure  20  we  plot  the  temperature  profiles  through  the  middle 
of  each  slice  (x  —  0)  at  the  final  time  t  —  3.46s  for  two  different  silicone 
thermal  conductivities.  The  top  plot  displays  temperature  profiles  for  the 
low  silicone  conductivity  ( ks  —  0.02 W/mK)-  the  bottom  plot  displays  tem¬ 
perature  profiles  for  the  high  silicone  conductivity  ( ks  —  1.02 W/mK).  Note 
that  with  a  low  silicone  conductivity,  the  larger  volume  fractions  result  in 
lower  temperatures  through  the  middle  of  the  slice.  The  total  temperature 
difference  between  profiles,  from  25%  to  55%  volume  fraction  at  the  bound¬ 
ary  74,  is  approximately  50 K.  The  reverse  behavior  is  true  in  the  case  of  a 
high  silicone  thermal  conductivity:  larger  volume  fractions  result  in  higher 
temperatures  through  the  middle  of  the  slice.  The  temperature  difference 
at  74  is  only  10 K. 

Again,  we  see  that  with  a  high  silicone  thermal  conductivity,  the  silicone 
slice  appears  to  be  more  efficient  extracting  the  heat  from  the  slice.  We 
refer  to  this  phenomena  as  the  speed  bump  behavior,  i.e.,  the  particles  in 
the  composite  slice  act  as  speed  bumps  and  slow  down  heat  transfer.  In 
Figure  21  we  plot  the  average  temperature  for  each  composite  slice  and  the 
silicone  slice  at  the  final  time  t  —  3.46s  across  a  range  of  silicone  thermal 
conductivities.  With  a  large  silicone  thermal  conductivity  ( ks  >  0.1  W/mK), 
larger  volume  fractions  result  in  higher  average  temperatures. 
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Average  temperatures:  t=3.46s,  kp=217W/mK 


Figure  22:  Average  temperatures,  comparing  heat  capacity 

3.3  The  Speed  Bump  Behavior 

The  speed  bump  behavior  seems  to  indicate  that,  at  the  higher  silicone  con¬ 
ductivity  range,  it  is  more  efficient  to  not  add  particles.  In  other  words,  the 
highly  conductive  particles  hinder,  rather  than  aid,  heat  dispersion.  This, 
of  course,  is  not  the  case.  The  speed  bump  behavior  is  due  to  the  different 
volumetric  heat  capacities  of  each  slice. 

The  all  silicone  slice  has  lower  temperatures  at  the  final  time  due  to 
the  fact  that  the  composite  slice  has  a  larger  volumetric  heat  capacity 
(pcp).  Though  starting  out  at  the  same  temperature,  the  slice  with  par¬ 
ticles  stores  more  heat  (energy)  than  the  all  silicone  slice.  In  the  previous 
pulsed  and  heated  slice  experiments,  the  aluminum  particles  have  a  heat 
capacity  of  pcp  =2.43xl06  J/m?K  whereas  the  silicone  has  a  heat  capacity 
of  pcp  =1.55xl06  J/m3K.  If  the  aluminum  particles  had  the  same  heat  ca¬ 
pacity  as  the  silicone  (i.e.,  same  stored  energy),  then  the  composite  slice 
with  highly  conductive  particles  should  be  as  or  more  efficient  at  dissipating 
heat  than  the  silicone  slice. 

To  test  this  we  adjusted  the  heat  capacity  of  the  aluminum  particles. 
In  Figure  22  we  plot  the  average  temperatures  at  the  final  time  t  —  3.46s 
across  a  range  of  silicone  thermal  conductivities  for  a  silicone  slice  (SIL),  and 
two  composite  slices  each  with  a  uniform  aluminum  particle  arrangement  of 
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25%  volume  fraction.  One  composite  slice  (UNI  hc-AL)  has  the  particle 
heat  capacity  of  aluminum  (i.e. ,  p  —  2.7e  —  12 g//r3,  cp  —  0.90 J/gK),  the 
other  composite  slice  (UNI  hc-SIL)  has  the  particle  heat  capacity  of  silicone 
(i.e.,  p  —  le  —  12 g/p37  cv  —  1.55 J/gK).  As  expected,  the  composite  slice 
containing  particles  with  the  heat  capacity  of  silicone  has  lower  average 
temperatures  than  the  all  silicone  slice.  In  other  words,  the  presence  of 
highly  conductive  particles  does  indeed  aid  the  heat  transfer. 

4  Continuous  Heat  Source  Experiments 

In  this  section  we  again  consider  the  basic  model  presented  in  Section  1, 
however  we  use  a  continuous  heat  source  at  the  boundary  74  rather  than  a 
laser  pulse  heat  source.  The  heat  source  is  a  constant  flux  term  of  the  form 
So(t )  =  Sc-  The  sides  71  and  73  are  insulated,  and  72  has  a  Newton  cooling 
boundary  condition.  See  Figure  2  for  the  pulsed  heat  source  experimental 
setup. 

For  completeness  and  comparison  purposes,  we  again  investigate  the 
following  effects  on  heat  transfer:  varying  particle  size  and  placement  for  a 
fixed  volume  fraction,  and  varying  volume  fraction  for  a  fixed  particle  size. 
To  this  end,  we  consider  the  steady  state  solutions  of  system  (1).  This  is 
obtained,  for  all  runs,  after  a  time  integration  of  100  seconds. 

4.1  Fixed  Volume  Fraction 

We  fix  the  volume  fraction  and  vary  the  size  and  placement  of  the  aluminum 
particles.  See  Figure  11  for  the  six  different  uniform  particle  arrangements. 

In  Figure  23  we  plot  the  temperature  profiles  through  the  middle  of  each 
slice  (x  —  0)  at  the  time  t  —  100s  for  two  different  silicone  thermal  con¬ 
ductivities.  The  top  plot  displays  temperature  profiles  for  a  low  silicone 
conductivity  (ks  —  0.02 W/mK)  slice;  the  bottom  plot  displays  temperature 
profiles  for  a  high  silicone  conductivity  ( ks  —  1  ,D2W/mK)  slice.  The  tem¬ 
perature  differences  are  slight,  still  the  composite  slices  with  smaller,  closer 
particles  obtain  lower  steady  state  temperatures  through  the  slice  middle. 
Also,  the  slices  with  particles  in  one  row  versus  two  rows  tend  to  channel  the 
heat  more  efficiently,  resulting  in  slightly  lower  steady  state  temperatures 
through  the  slice  middle. 

The  same  trends  are  observed  with  the  average  temperatures  for  various 
silicone  conductivities.  In  Figure  24  we  plot  the  average  temperature  for  each 
of  the  uniform  particle  arrangements  and  the  corresponding  silicone  slice  at 
the  final  time  t  —  100s  across  a  range  of  silicone  thermal  conductivities.  Note 


25 


t=1 00s,  k  =0.02  W/mK 

s 


y 


Figure  23:  Temperature  profiles  at  x  —  0,  25%  volume  fraction 


Figure  24:  Average  temperatures  for  slices  with  a  continuous  heat  source 


that  the  range  of  average  temperatures  is  larger  than  for  the  case  of  a  pulsed 
heat  source  (compare  with  Figure  19).  With  the  lowest  silicone  conductivity, 
average  temperature  differences  between  silicone  and  composite  slices  are  as 
much  as  20 K.  With  a  high  silicone  conductivity,  the  average  temperature 
difference  drops  to  IK. 

4.2  Varied  Volume  Fraction 

In  this  section  we  look  at  the  effect  of  volume  fraction  on  heat  transfer 
through  a  composite  adhesive  with  a  continuous  heat  source.  The  particle 
size  is  fixed  ( R  —  10.50/i),  the  volume  fraction  varies  from  25%  to  55%. 
The  particles  are  aluminum  (kp  —  217W/mK)  and  arranged  randomly.  See 
Figure  14  for  the  random  particle  arrangements. 

In  Figure  25  we  plot  the  temperature  profiles  through  the  middle  of  each 
slice  ( x  —  0)  at  the  final  time  t  —  100s  for  two  different  silicone  thermal 
conductivities.  The  top  plot  displays  temperature  profiles  for  a  low  silicone 
conductivity  ( ks  —  0.02 W/mK)  slice;  the  bottom  plot  displays  temperature 
profiles  for  the  high  silicone  conductivity  ( ks  —  1.02  W/mK)  slice.  With 
both  conductivities,  the  slices  with  larger  volume  fractions  obtain  lower 
temperatures  through  the  middle  of  the  slice.  We  obtain  similar  results  as 
the  pulsed  heat  source  scenario;  compare  with  Figure  20. 
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Finally,  in  Figure  26  we  plot  the  steady  state  average  temperatures  for 
each  of  the  volume  fraction  slices,  and  the  all  silicone  slice,  at  the  final  time 
t  —  100s  across  a  range  of  silicone  thermal  conductivities.  The  trend  is  the 
same  as  in  the  pulsed  heat  source  scenario:  the  larger  the  volume  fraction 
the  lower  the  average  temperature.  Additionally,  all  the  composite  slices 
obtain  a  lower  average  temperature  than  the  all  silicone  slice  for  the  full 
range  of  silicone  thermal  conductivities. 

5  Summary 

We  investigated  the  heat  transfer  through  composite  adhesives  under  three 
experimental  scenarios:  a  laser  pulse  heat  source,  a  preheated  sample,  and  a 
continuous  heat  source.  In  the  pulse  heat  source  scenario,  we  reported  on  a 
variety  of  results  (average  flux,  average  sensitivity)  to  better  understand  the 
basic  heat  transfer  behavior.  Similar  results  are  available  for  the  preheated 
and  continuous  heat  source  scenarios  as  well. 

For  all  three  scenarios  we  considered  the  effect  of  particle  size  and  place¬ 
ment  for  a  fixed  volume  fraction  as  well  as  the  effect  of  volume  fraction 
for  a  fixed  particle  size.  Results  for  the  continuous  and  pulse  heat  sources 
were  similar.  Namely,  smaller,  closer  particles  are  slightly  more  effective  at 
transferring  heat  through  a  composite  slice.  Composites  with  higher  volume 
fractions  are  also  more  effective.  The  effectiveness,  however,  depends  on  the 
thermal  conductivity  of  the  silicone.  In  other  words,  with  a  large  silicone 
thermal  conductivity,  these  effects  become  less  pronounced.  Alternatively, 
when  working  with  composites  with  low  silicone  thermal  conductivities,  the 
heat  transfer  may  be  improved  by  slightly  increasing  the  silicone  thermal 
conductivity. 

These  conclusions  apply  to  the  preheated  scenario  as  well,  but  are  not 
evident  in  the  temperature  profiles  and  average  temperature  results  dis¬ 
played  in  Section  3.  We  informally  called  this  discrepancy  the  speed  bump 
behavior.  It  is  due  to  the  fact  that  the  preheated  composite  slices  have  a 
larger  volumetric  heat  capacity  than  an  all  silicone  slice.  In  other  words, 
more  heat  is  stored  and  transfered  in  the  composite  slices.  Similarly,  larger 
volume  fraction  slices  have  a  larger  volumetric  heat  capacity. 
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